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Abstract 

We explore the impact of the data from the KamLAND experiment in constraining neu- 
trino mass and mixing angles involved in solar neutrino oscillations. In particular we discuss 
the precision with which we can determine the the mass squared difference Attt-q and the 
mixing angle ^0 from combined solar and KamLAND data. We show that the precision 
with which Atuq can be determined improves drastically with the KamLAND data but the 
sensitivity of KamLAND to the mixing angle is not as good. We study the effect of enhanced 
statistics in KamLAND as well as reduced systematics in improving the precision. We also 
show the effect of the SNO salt data in improving the precision. Finally we discuss how a 
dedicated reactor experiment with a baseline of 70 km can improve the 9q sensitivity by a 
large amount. 
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1 Introduction 



Two very important results in the field of neutrino oscillations were declared in year 2002. The 
first data from the neutral current (NC) events from Sudbury Neutrino Observatory(SNO) exper- 
iment were announced in April 2002 [1] Comparison of the the NC event rates with the charged 
current (CC) event rates established the presence of v^/vr component in the solar flux reinforc- 
ing the fact that neutrino oscillation is responsible for the solar neutrino shortfall observed in the 
Homestake, SAGE, GALLEX/GNO, Kamiokande and SuperKamiokande experiments. The global 
analysis of solar neutrino data picked up the Large Mixing Angle (LMA) MSW as the preferred 
solution [2] . The smoking-gun evidence came in December 2002 when the KamLAND experiment 
reported a distortion in the reactor anti-neutrino spectrum corresponding to the LMA parameters 
[3] . The induction of the KamLAND data in the global oscillation analysis resulted in sphtting 
the allowed LMA zone in two parts (at 99% C.L.) - low-LMA lying around Am| = 7.2 x 10~^ eV^, 
sin^ ^0 = 0.3, and high-LMA with AttIq = L5 x lO^"' eV^, sin^ = 0.3 respectively. The low- 
LMA solution was preferred statistically by the data [4]. The recently announced SNO data from 
the salt-phase [5] has further disfavoured high-LMA and it now appears at > 99.13% C.L. [6]. 
Thus the SNO and KamLAND results have heralded the birth of the precision era in the measure- 
ment of solar neutrino oscillation parameter. In this article wc take a closer look at the precision 
with which we know the solar neutrino oscillation parameters at present and critically examine 
how precisely they can be measured with future data. 

2 Oscillation Parameters from solar neutrino data 

In fig. 1 we show the impact of the SNO NC data from the pure D2O phase, the salt phase as well 
as combining the information from both phases on the oscillation parameters Ato|, (= Am|i ) and 
sin^ 9q{= sin^ 612) from a two-flavour analysis . We include the total rates from the radiochemical 
experiments CI and Ga (Gallex, SAGE and GNO combined) [7] and the 1496 day 44 bin SK 
Zenith angle spectrum data [8]. For the pure D2O phase we use the CC-f-ES-f-NC spectrum data 
whereas for the salt phase we use the published CC,ES and NC rates [6]. The details of the 
analysis procedure can be found in [9]. Also shown superposed on these curves are the isorates of 
the CC/NC ratio. We find that 

• The upper limit on Am|, tightens with the increased statistics when the salt data is added to 
the data from the pure D2O phase. 

• The upper limit on sin^ tightens. For the neutrinos undergoing adiabatic MSW transition 
in the sun Rcc/Rnc ~ sin^ 9. The SNO salt data corresponds to a lower value of the CC/NC 
ratio which results in a shift of sin^ ^0 towards smaller values. 

3 Impact of KamLAND data on oscillation parameters 

The KamLAND detector measures the reactor antineutrino spectrum from Japanese commercial 
nuclear reactors situated at a distance of 80 -800 km. In this section we present our results of 
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Figure 1: The 90%, 95%, 99% and 99.73% CL. allowed regions in the Am2i — sin^ 6q plane from 
global x^-analysis of the data from solar neutrino experiments. We use Ax^ values to plot the 
CL. contours corresponding to a two parameter fit. Also shown are the lines of constant CC/NC 
event rate ratio Rcc/nc- 

global two-generation analysis of solar+KamLAND spectrum data. For details we refer to our 
analysis in [4, 10] . 

Figure 2 shows the allowed regions obtained from global solar and 162 Ton-year KamLAND 
spectrum data. As is seen from the leftmost panel of figure 2 the inclusion of the KamLAND data 
breaks the allowed LMA region into two parts at 99% C.L.. The low-LMA region is centered 
around a best-fit Am^ of 7.2 x 10~^ eV^ and the high-LMA region is centered around 1.5 x 10~^ 
eV^. At 3cr the two regions merge. The low-LMA region is statistically preferred over the high- 
LMA region. With the addition of the SNO salt data the high-LMA solution gets disfavoured 
at 99.13% C.L.. In Table 2 we show the allowed ranges of Am^ and sin^ 6'© from solar and 
combined solar+KamLAND analysis. We find that ^rn^^ is further constrained with the addition 
of the KamLAND data but sin^ 612 is nor constrained any further. 
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Tabic 1: The best-fit values of the solar neutrino oscillation parameters, obtained using different 
combinations of data sets. Shown also are the 99% G.L. allowed ranges of the parameters from 
the different analyses. 
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Table 2: The range of parameter values allowed at 99% G.L. and the corresponding spread. 



4 Closer look at KamLAND sensitivity 

In Table 4 we take a closer look at the sensitivity of the KamLAND experiment to the parameters 
Am2i and 612 with the current as well as simulated future data and examine how far the sensitivity 
can improve with the future data. We define the % spread in oscillation parameters as 

spread = x 100 (1) 

and determine this quantity for the current solar and KamLAND data as well as increasing the 
KamLAND statistics. The current systematic error in KamLAND is 6.42% and the largest contri- 
bution comes from the uncertainty in fiducial volume. This is expected to improve with the calibra- 
tion of the fiducial volume and we use a 5% systematic error for 1 kTy simulated KamLAND data 
and 3% systematic error for 3 kTy simulated KamLAND data. The table reveals the tremendous 
sensitivity of KamLAND to Amli- The addition of the present KamLAND data improves the 
spread in Am2i to 30% from 68% obtained with only solar data. With 1 kTy KamLAND data it 
improves to 10% and if we increase the statistics to 3 kTy then the uncertainty in Amli reduces 
to 6%. However the sensitivity of KamLAND to the parameter 612 does not look as good. The 



4 



lO 



' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

CI -1- Ga. -+■ Sk spec 
-+- KamLAND 



lO 




I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 

Cl H- Ga ■+- SK spec 
-H SNO Csalt) 
-+- KamLAND 




I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 

Cl -H Gra -H Sk spec 

H- SNO Csalt) 
ICamLAND 




5 I, , I , , , , I , , , , I , , , , I , ,1, , I , , , , I , , , , I , , , , I , ,1, , I , , , , I , , , , I , , , , I , 

0.2 0.3 0.4 0.5 0.2 0.3 0.4 0.5 0.2 0.3 0.4 0.5 



- 2 

sm fc) 



1 2 



Figure 2: The 90%, 95%, 99% and 99.73% CL. allowed regions in the Amgi — sin^ 9q plane from 
global x^-analysis of solar and KamLAND data. We use the Ax^ values corresponding to a 2 
parameter fit to plot the CL. contours. 

addition of the current KamLAND data to the global solar analysis does not improve the spread in 
sin^ 9i2- With reduction of the systematic error to 5% the spread with 1 kTy statistics improves to 
26% and even with a very optimistic value of 3% for the systematic uncertainty and a substantial 
increase of statistics to 3 kTy, the KamLAND data fails to constrain 612 much better than the 
current solar neutrino experiments. 

In figure 3 we compare the allowed areas computed with spectrum simulated at Am^i = 
7.2 X 10~^ eV^ and Am^i = 3.5 x 10~^ eV^. We show limits for the current KamLAND systematic 
uncertainty of 6.42% and a very optimistic systematic uncertainty of just 2%. The % spread in 
uncertainty for the spectrum simulated at 7.2 x 10~^ eV^ with 6.42% systematic uncertainty is 
37% while for 3.5 x 10^^ eV^ case the spread is 25%. The effect of reducing the systematics to 2% 
results in the spread coming down to 32% and 19% respectively [11]. We would like to mention 
that the figure 3 uses the CC, NC and ES rates from the D2O and not the latest results from the 
salt phase. However the purpose of this figure is to compare the spread in sin^ 6q obtained for the 
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Figure 3: The contours for the combined analysis using the solar and 3 kTy simulated KamLAND 
spectrum. The first two rows of panels correspond to spectrum simulated at 7.2 x 10~^ eV^ while 
the lowermost row of panels are for KamLAND data simulated at Amli = 3.5 x 10^^ eV^. 
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Figure 4: The probability vs distance for an average energy of 4 MeV for Am^i = 7.5 x 10 ^eV^ 
and 3.5 xlQ-^ eV^. 

two different values of ^m\^ and the use of the sno salt phase data is not going to change the 
relative spreads significantly. We also present in the middle panel of 3 the allowed areas drawn 
using a 7% uncertainty in the NC rate. The uncertainty in the NC rate from the 1^20 (salt) phase 
data is 12% (9%). 

To trace the reason why the 9 sensitivity is better at 3.5 x 10^^ eV^ in figure 4 we plot the 
probability vs distance for energy fixed at 4 MeV. The figure shows that the average distance of 
~ 150 km of KamLAND corresponds to a maximum in the probability for 7.2 x 10~^ eV^ while 
at 3.5 X 10~^ eV^ corresponds to a minima. 

The relevant survival probability for KamLAND is given by the vacuum oscillation expression 




(2) 



where Li stands for the different reactor distances and one needs to do an averaging over these. 
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Figure 5: The survival probability Pee as a function sin^6'i2. 
Three limits can be distinguished 

• sin^(Am^i L/AE) = we get a Survival Probability MAXimum (SPMAX) 

• sin^(Am^i L/AE) = 1 we get a Survival Probability MINimum (SPMIN) 

• sin^(Am2i L/AE) = 1/2 we get averaged oscillation 

In the LMA region the solar neutrinos undergo adiabatic MSW transition and the survival 
probability can be approximated as 

Pee('i? )~sin2^0 (3) 

Whereas the low energy pp neutrinos do not undergo any MSW resonance and the survival prob- 
ability is just the averaged oscillation probability in vacuum. 

In figure 5 we plot the 9 dependence of the adiabatic MSW probability as well as the probability 
for the SPMIN and averaged oscillation case [11]. The figure shows that for large mixing angles 
close to maximal, the adiabatic case has the maximum sensitivity. For mixing angles not too 
close to maximal (sin^ ^0 < 0.38) , the Pee for the SPMIN case has the sharpest dependence on 
the mixing angle and the 6q sensitivity is maximum. Since the 99% CL. allowed values of 6q is 
within the range 0.22 < sin^ ^0 < 0.44, SPMIN seems most promising for constraining Ou- On 
the other hand at SPMAX the oscillatory term goes to zero and the 612 sensitivity gets smothered. 
Since in the statistically significant region the KamLAND probability corresponds to an SPMAX 
for he best-fit value of 7.2 x 10~^ eV^ the 6 sensitivity of KamLAND is not as good as its Arn^i 
sensitivity. For this value of Am^i the SPMIN comes at 70 km. 



8 



4 5 e> 7 



4 5 e. 7 « 

E^^3 Cin MeV) 



lO 



— 1 1 1 — I — I — I I I I 1 1 1 — I — I — I I 1 

NEW : 



-9^ 

lO" 



lO 



-9^ 

-HcvT:: lo" 



— 1 1 1 — I — I — I I I I 1 1 1 — I — I — I I 1 

Solar -+- NEW : 



lO 



_l I I I I I I 



_l I I I I I I 



lO 



_l I I I I I I 



_l I I I I I I 



lO 



lO 



lO 



lO 



lO 
tan^e 



lO 



Figure 6: The spectrum and allowed regions for a new reactor experiment with a source detector 
distance of 70 km 

5 A dedicated reactor experiment for sin^ ^0 

We show in Figure 6 the constraints on the mass and mixing parameters obtained using a "new" 
dedicated reactor experiment whose baseline is tuned to an oscillation SPMIN [11]. We use the 
antineutrino flux from a reactor a la Kashiwazaki nuclear reactor in Japan with a power of about 
25 GWatt. We assume a 80% efficiency for the reactor output and simulate the 3 kTy data at 
the low-LMA best-fit for a KamLAND like detector placed at 70 km from the reactor source and 
which has systematic errors of only 2%. The top-left panel of the Figure 6 shows the simulated 
spectrum data. The histogram shows the expected spectrum for no oscillations. E^s is the 
"visible" energy of the scattered electrons. The top-right panel gives the ratio of the simulated 
oscillations to the no oscillation numbers. The sharp minima around 3 — 4 MeV is clearly visible. 
The bottom-left panel gives the CL. allowed areas obtained from this new reactor experiment 
data alone. With 3 kTy statistics we find a marked improvement in the 6q bound with the 99% 
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Figure 7: The isorates for a pp scattering experiment. 



range 0.39 < tan^ 6q < 0.52 giving a spread of 14% . The "dark side" solution appearing in the 
left lower panel because of the 9q — {-k /2 — 6q ) ambiguity in the vacuum oscillations probability 
is ruled out in the right lower panel by the solar neutrino data. Recently sites of reactor neutrino 
experiments with a source- reactor distance of 70 km has been discussed in [12]. Also in Japan a 
new reactor complex SHIKA-2 at ~ 88 km (close to SPMIN) will start in 2006 (See however [13]). 



6 Other future experiments 

In Figure 7 we show the lines of constant rate/SSM predicted in a generic LowNu electron scat- 
tering experiment sensitive to pp neutrinos [14]. At these low energies neutrinos the survival 
probability in the LMA zone is ~ 1 — | sin^ 2^ and has almost no sensitivity to /S.m'^^. But the 
sin^ Oq sensitivity is quite good and thus these experiments may have a fair chance to pin down 
the value of the mixing angle ^12, if they can keep the errors low. 
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7 Conclusions 



The KamLAND experiment has not only confirmed the LMA solution to the solar neutrino prob- 
lem it has narrowed down the allowed range of Am| considerably owing to its sensitivity to the 
spectral distortion driven by this parameter. However the 6q sensitivity of KamLAND is not as 
good. Baseline is important to identify which parameters would be best determined. We discuss 
that a SPMIN in the vacuum oscillation probability is important for the determination of the 
mixing angle. For the current best-fit Amg in the low-LMA region SPMIN comes at a distance 
of 70 km We propose a dedicated 70 km baseline reactor experiment to measure 9q down to 
~ 10% accuracy. LowNU experiments could be important for precise determination of 6q if the 
experimental errors are low. 

This talk is based on the work [11]. The updated analysis including the SNO salt results were done 
in collaboration with S.T.Petcov and D.P. Roy and the authors would like to acknowledge them. 
This work was supported in part by the Italian MIUR and INFN under 'Fisica Astroparticellare" 
(S.C.). 
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